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Using many-body theory we develop a set of formally exact kinetic equations for inhomo- 
geneous condensate and one-body observables. The method is illutrated for (p^ field theory in 



^ , 1+1 dimensions. These equations, when computed with the help of time-dependent projec- 

Zj . tion technique, lead to a systematic mean-field expansion. The lowest and the higher order 

terms correspond to, respectively, the gaussian approximation and the dynamical correlation 

/\ ' effect. 
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I. Introduction 

The dynamical evolution of inhomogeneous field configurations is an important problem 
and a common theme in cosmology, high energy and condensate matter physics. In cos- 
mology inhomogenous field configurations appear when topological objets such as textures 
or cosmic strings involving inhomogenous field configurations are considered. Their relax- 
ational dynamics is thought to have a bearing on the spectrum of fluctuations in the cosmic 
microwave background radiation [|I[]. In the ultrahigh energy heavy-ion collisions {^/s > 200 
GeV/nucleon) a large energy density (few GeV/fm^) is deposited in the collision region cor- 
responding to temperature above the critical value for chiral symmetry restoration. In this 
situation, it is possible within a volume of few fm^ an inhomogeneous condensate is formed. 
When these regions cool down this field relax toward the equilibrium situation which might 
be misaligned with the vacuum state 0. On the other hand, the recent sucess of the experi- 
mental observation of Bose- Einstein condensation for systems of spin polarized magnetically 
trapped alkah atoms at ultra-low tempeture is smulating development of theory for the evo- 
lution of nonuniform condensates . 

A microscopic description of this off-equilibrium process require a nonperturbative treat- 
ment of inhomogeneous condensates and in the field-theoretical context, this has been im- 
plemented through the use of a Gaussian ansatz for the wavefuncctional in the framework of 
a time-dependent variational principle @]. Actually, the Gaussian ansatz, having the form of 
an exponential of quadratic form in the field operators, implies the many-point functions, can 
be infact factord in terms of two-point functions. The dynamics of the reduced two-point 
density becomes then itself isoentropic, as a result of irreducible higher-order correlation 
effects being neglected. 

The purpose of the present paper is twofold. The first one is to reevaluate and improve 
the gaussian approximation. We follow a time-dependent approach developed earlier in the 
context of nuclear many-body dynamics 0. This methos allows for a formulation of a 
mean-field expansion for the dynamics of the two-point correlation function from which one 
recover the results of the gaussian mean-field approximation in lowest order. Beyond this, 
we are able to explicitly incluide higher dynamical correlation effects. The second purpose is 
to extend the former results in the context of spatial uniformity to the inhomogeneous field 
configurations P] . In this case the spatial dependence of the field operator is expanded in the 
general natural orbitals. These orbitals can be given in terms of an expansion of convenient 
basis which willm also evolve in time according to additional dynamical equations. Although 
the procedure is quite general, we will apply our method in the context of a single scalar 



field in 1+1 dimensions. This will illutrate all the relevant points of the approach and cut 
down inessential technical complications. 

II. A Simple Example in Quantum Mechanics 

This section will illustre the scheme of approximation in the context of 0^ in + 1 
dimension 0. Although this theory is finite, we introduce nevertheless a counterterm for 
the further use. Therefore the hamiltonian reads as 

2 2 ^ 24^ 8m^ ^ ^ 

where m is the renormalized mass. In the usual quantum mechanics (QM) language, one 
speaks of a particle in the quartic potential with (j) and n being its position and momentum 
operator and satisfy the usual communtation relation, [0, vr] = i. In the quantum field theory 
(QFT) language, however, one imagines that this field lives in one point and the particles 
have quartic self-interaction [Q. 

A natural choice of the subsistem in the context of many-body problem is the observables 
associated to the one-body density. The exact microspic description of the time evolution of 
the one-body density of a many-body system can be formally given in terms of a sum of two 
parts: the usual gaussian contribution (also known as time- dependent hartree-bogoliubov 
approximation) plus additional dynamics arise from the time evolution of quantum correla- 
tions in the entire system. The physical origin of the later contribution lies in the complicated 
dynamical evolution of quantum correlationin the entire system and their consequences in 
the change in the coherence properties of each system. 

Il-a. Gaussian Variables 

In order to derive the dynamical equations for the one-body observables we focus on the 
operators which are either linear or bilinear forms of creation and annihilation, henceforth 
referred to as gaussian observables. We begin therefore write the Heisenberg field operators 

and as 

0(t) = -^ [a\t) + a(t)] n{t) = i.l^[a\t) - a(t)] (2.2) 

where a, a) are the usual annihilation and creation operators satisfying the boson commuta- 
tion relations: [a, a^] = 1; the parameter jj, will be fixed later in a convenient way. The state 
of the entire system is given in terms of the matrix density F in the Heisenberg picture. It 
is Hermitean, time indenpedent and has a unit trace. 



The first variable of interest is given as mean value of annihilation, 

At = Tra(t)F, (2.3) 

and its complex conjugate. From this one can define the shifted boson operator 

b = a-A. (2.4) 

Next, we consider the possible mean values of the bilinear forms of the shifted operators, 
which can be combined in an extended one boson plus pairing density, 

\ (6t(t)6t(t)) {h{t)h\t)) I \W i + A*y 

The quantities A and 11 together with the mean-value of field. A, describe the one-body 
observables of the system. 

To deal with the pairing density 11 we proceed, as usual, by defining the Bogohubov 
quasi-particle operator as 

d{t) = xlh{t) + ylh\t) S{t) = xtb^it) + vMt) (2.6) 

and require that (dd) = {d'^d'^) = 0. A sistematic way to determine the coeffcients of the 
Bogoliubov transformation, xt and yt, is to solve the foUwoing secular problem: 

GRX = XGN , (2.7) 

where 

O^i' M . xJ^'A . A^ J"' » 1 . ,2.8) 

V -1 y \yt X* I \^ 1 + i^i y 

The eigenvalues vt = Trd^{t)d{t)F stand for occupation numbers of quasi-particles. Since 
the Bogoliubov transformation is canonica one can verify that X^ satisfies the orthogonality 
and completeness relation, i.e., 

X^GX = XGX^ = G (2.9) 

In terms of coeficients of transformation (2.9) reads as 

\xt\^-\yt |'= 1 (2.10) 



Thus, the paring densities {d{t)d{t)) and {S{t)S(t)) can be obtained from xt, yt and ut. The 
next step is to get equations of motion for these quantities. 

Il-b. Equation of Motion and Gaussian Approximation 

The Next step isto obtain the time evolution for the variables of interest discussed in 
the previous subsection. We begin with the amplitude of condensate At defined by Eq.(2.3). 
Using Heisenberg equation of motion one has 

iAt = Tr[a, H]F = XiTr[d, H]F + Tvy*[d\ H]F (2.11) 

where we have used (2.4) and (2.6). Now, to get equantions for the Bogoliubov coeficients 
we rewrite (2.7), using (2.9), as 

XJRXt = N. (2.12) 

Taking time derivative of this and using Eq.(2.9) we get 

XJRXt = N - XJRXt + Xj RXt ■ (2.13) 

The left-hand side of this equation can be evoluated from the Heisenberg equation of motion 
and (2.6) and yields 



iXlRXt 



Tr[d^d,H]F TT[dd,H]F 



±1IU U..JJU' ±L\UU.1I\1' \ 

(2.14) 
\ TT[d^d^F Tr[dd\H]F 



The right-hand side of (2.13), on the other hand, can be obtained with the help of (2.9). 
Equation now the result to (2.14) one gets 

iiy = TT[dU,H]F (2.15) 

and 

i{xy - xy) = Tr [dU, H]F. (2.16) 

The equations (2.10), (2.15) and (2.16), together with the normalization condition, deter- 
mine fully, in principle, the time evolution for the gaussian observables if the density matrix 
F is expressible in terms of the quantities themselves. However, when the hamiltoniana H 
involves self-interaction fields, traces of these equations will involve also many-body densities 
and therefore they are not closed. One emergent approximation to deal with this situation 
is to replace the full density F by a truncated one, Fq, which has the form of a exponential 



of bilinear in creation and annhilation operators (see, e.g., Eq.(2.22) of 0). In terms of 
quasi-particular operator Fq can be written conveniently as 

1 / Z/, \dHt)d{t) 

Fo{t) = — — -^) (2.17) 

Notice that Fq is diagonal in the quasi-particle basis and contain no irreducible two or many 
quasi-particle correlations. Furthemore, it is easy to verify that Fq given by (2.17) has a 
unit trace and reproduce the corresponding average of linear and bilinear field oprarors of 
full density, i.e., 

Tr aFQ = A = Tr aF 
Tr SdFQ = // = Tr SdF 
Tr ddFQ = = Tr ddF 

Therefore the approximation gives a set of self-consistent equations for the one-body observ- 
ables and is is what we refer as the mean-field or gaussian approximation. In particular, when 
Fq is written in terms of field representation, it is equivalent to the density used by Jackiw 
@]. This approximation constrains the time evolution of the system to remain in a gaussian, 
which contains no irreducible two or many particle correlation, whereas the true evolution 
will, as time progresses, introduce (quasi)particle-particle correlation not describable by the 
gaussian-like matrix density. In terms of general discussion made in the previous subsection 
the limitation of the approximation shows up in the dynamical evolution of occupancy. 



ZP 



Tr [dU, H]Fq = Tr [dU, Fq]H = 0. (2.18) 



Therefore , further improvements have to be achieved in order to describe correlation effects 
between different subsistems. 

II-c. Projection Technique and Dynamical Correlations 

The question we want to address ourselves at this point is how to express the full cor- 
related density F of entire system in terms ingredients of subsystem (gaussian variables) in 
such way that one can get a set of closed equations for the one-body observable. A framework 
to achieve this goal was developed some time ago by Willis and Picard using time-dependent 
projection technique |jlO| in the context of master equation for coupled systems. The method 



was extented later by Nemes and Piza to study nuclear many-body dynamics f^ . The method 
consists essentially in writing the correlation informations of the full density in teerms of a 
memory kernel acting on the uncorrelated density Fq . 
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Following their strategy we first decompose F in two parts 

F = Fo{t) + F'{t) 



(2.19) 



where -fo(^) is the uncorrelated part given by (2.17) and therefore F'{t) is a traceless 
correlation part. The substitution of F by just Fo(t) in the equations (2.11) and (2.15)- 
(2.16) gives the usual gaussian approximation as we have discussed before. The crucial step 
is to observe that Fo(t) can be seen as a time- dependent projection of F , i.e., 

Fo = P(t)F , P(t)P(t) = P(t) . (2.20) 

For the explicit construction of P(t) we require, in addition to eqs. (2.20), the condition 

zFoit) = [Fo{t),H] + F{t)[H,F] (2.21) 

which is the Heisenberg picture counterpart of the Schrodinger picture condition used in ||n 



to determine P uniquely. The resulting form for P(t) is (see appendx A of Ref. for 
details of the derivation) 



P- 



+ 



1- 
dd 



dU 



P 



dU 



d 



S 



l+P 
Tr (S S •) 



Tr(-) + " \ Trfrftrf.) + _ Tr f rft .) + -^1— Tr(rf 

p(l + p) ^ ' [p ^^ l-|-p 



Ti (dd- 



(2.22) 



2pp ' V~" " / ■ 2{l+p){l+p) 

where the dot stands for objets on which the projector acts. 

The next step is to obtain a differential equation for the correlated density F'{t). This 
follows immediately from eqs. (2.19) and (2.21), 

/. d 



I- -P(t)P F'(t) 



Q(t)PFo(t) , 

where we have introduced the operators 

Q(t) = IL-P(t) , P- = [iJ, ■] . 

This equation has the formal solution 

F'it) = G(t,0)F'(0) - i f dt'G(t,t')m')I.Fo(t'] 

Jo 

where G{t, t') is the time-ordered Green's Function 

G{t,t') = T expi f dT'P{T)L . 
Jf 
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(2.23) 



(2.24) 



(2.25) 



(2.26) 



What we have obtained so far is a formaly exact expression relating F'{t) and Ftr (for t' <) 
and the initial correlations F'{0). This allows us to get a set of closed dynamical equations 
as traces over functional of -fo(^') and initial correlations. An actual calculation is, however, 
hopeless because of the complicated time dependence of the Heisenberg operators present in 
the memory kernel of (2.25). A systematic expansion to treat the memory integral of the 
equation (3.8) has been discussed in Ref . [pJj] in the Schrodinger picture. The implementation 
of the corresponding expansion in the Heisenberg picture consists in approximating the time 
evolution of the field operators, when evoluating memory effects, by a simpler one-body 
generator of mean-field evolution, i.e.. 



id = [d, Hq] — iA + i{xlxt — illyt)d 



■ / • * * 



yWd^- 



(2.27) 



The last three terms account for the (explicit) time dependence of d{t) on condensate and 
pairing effects; The hamiltonian Hq is taken, in this approximation, as a mean-field one, 

Hq = P^ H + d'^ Tr[d ,H]F' - dTi \S , h] F' 



rftrft 



2(1 + 2p) 



Tr[dd, H]F' 



dd 



2(1 + 2p) 



Tr 



d^ S , h] f'. . 



(2.28) 



This approximation results a unitary time evolution for the Heisenberg field operator and 
operators in different time are related by a phase factor 



d{t) 



jip{t,t') 



d{t') 



(2.29) 



(see the third equation on page 1609 of 0; see also (5.3) of |^ for details). 

In this way, the authors of the ref.|jll| devised a systematic expansion for the correlation 
density F'{t) [see their Eq.(3.10)]. The corresponding expansion in the Heisenberg picture 
reads as 



F'{t) = (G(t,0)F'(0) - I f <qti)lLFo(ti 

Jo 



dti 



dt2Q;t2)(lL-lLo(t2) 



Q;ti)lLFo(ti) + 



(2.30) 



where ILq ■ = [Hq , ■]. In what follows we restrict ourselves to initial conditions such that 
-F'(O) = and to the lowest approximation for F'{t). Therefore, the evaluation of the 



equations of motion involves traces of the type 

Ti[d{t),H]Fo{t) - iTi[d{t),H] f'dt'mt')[H,Fo{t')] 

Jo 

mean — field correlation (2.31) 

where 0{t) can be d{t) , d'^{t) d{t) , d{t) d{t) , and operators at different times are related 
by equation (2.29). 

We have now all the necessary ingredients for implementation of the approximation and 
the derivation for the equation of motion is a straightforward algebraic exercise. The results, 
including the numerical calculation, for (f)Q_^_i model are shown in ^. Extensions of this 
method to other models of field theory were obtained recently . The results demonstrate 
that this approach can overcome some conceptual dificults of gaussian approximation as well 
as a much better description for the gaussian observables. 

III. Kinetic Equations for Simple Observables of the Field 

Section II reviewed some important points of our extended gaussian approximation and 
its implementation in the simplest context of quantum mechanics. In this and next sec- 
tion we will report its applications in context of inhomogenoeus field configuration, which 
is relevant for the dynamical evolution of many body finite system. The simpler case of 
spatial uniformity was discussed for sevaral field models recently. Technical dificulties in 
such cases reduce tremendously because of translational invariance and the gaussian vari- 
ables are automatically diagonal in momentum space. In other words, the eigenfunctions 
of the body-density, known as natural orbitals, are independent of time and given by plane 
wave. For this new scenary, however, the natural orbitals are time dependent. Therefore, 
additional equation is needed in order to get self-consistent equations of motion. 

Ill-a. Generalized Bogoliubov Transformation 

In order to implement the idea let us first expand the Heisenberg field operator (j){t, x) 
and the canonical momentum 7i{t, x) as 

(P{t,x) = 5: [/fc(x)afc(t) + /;(x)4(t)] , (3.1) 

k 

7r(t,x) = -iJ2ko .fk{,x) akit) - fl{x) a\{t) , (3.2) 



where ak{t), al{t) are boson operators satisfying the equal time commutation relation 

[ak{t) , ak'{t)] = 5kk' ■ (3.3) 

The fk{.x) are the periodic boundary condition plane waves 



^-^tfZ ' lL 



L being the lenght of the periodicity box and k"^ = fc^ + //^ . The expansion mass 
parameter fi is conveniently fixed, e.g. in terms of the equilibrium solution in the mean 
field approximation p. 

The next step is to focus on the variables of interest, which are mean- value of linear and 
bilinear boson operators. The first of them is the expectation value of the field operator, 

{(pit,x)) = Ti(j){t,x) F , (3.5) 

where F is a density matrix in the Heisenberg picture that characterizes the state of the 
system. In terms of the expansion (2.1), one has 

(0(t,x)) = Y: ihi^) Mi) + fm Al{t)] (3.6) 

k 

with 

Akit) = Tiak{t)F . (3.7) 

We can now define the shifted boson operators with the help of the Akit), 

bk{t) = ak{t) - Ak{t) , (3.8) 

and include as variables of interest also the expectation value of pairs of bk{t), 6|.(t) at 
equal times: 

Akk'it) = Tibl,{t)bkit)F , (3.9) 

Ukk'it) = Tibk'{t)bkit)F . (3.10) 

The hermitean matrix A and the symmetric matrix 11 are in fact the one-boson density 
matrix and the pairing density for the shifted bosons respectively. The corresponding 
matrices for the a-boson are easily expressed in terms of A , 11 and Ak . 
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The next step is to write the one-body density matrix in diagonal form and also incor- 
porate information of the pair density in an associated set of natural orbitals by setting up 
an extended one-body density matrix 0, 

(a n \ 

R = \ \ = R^ (3.11) 

\n* i + A* J 

and solving the extendend version of eigenvalue problem defined by (2.7), where the matrix 
elements are given now the following matrices: 



1 \ U V* 

G = \ \ , X = \ \ , N 

-1 / \ V u* 



( P 

(3.12) 

V 1 + P 



Since (2.7) is a non-hermitean eigenvalue problem, it is important to consider also the adjoint 
equation 

RGX = XGN , (3.13) 

from which one finds that 

X = GX . (3.14) 

The adjoint vectors X satisfy biorthogonality relations with X which allow one to introduce 
the normalization condition 

X^X = X^GX = G (3.15) 

and the completeness relation 

XGX^ = G . (3.16) 

Furthermore, one can use the eigenvectors of the secular problem (2.12) to construct new 
boson operators 

^ d \ , ( b . 

Xt . (3.17) 



V 



This equation defines in fact the general Bogoliubov transformation 

da{t) = E (uUt) hit) + V:^{t) blit)) . (3.18) 

k 

It is easy to see that the secular problem (2.12) constrains the expectation value of products 
of da(t), dl^{t) as 

Tr 4(t) d^it) F = P^{t) 5afs (3.19) 

Tr d^t) df3{t) F = . (3.20) 

11 





where Pa are elements of the eigenvaue matrix P and can be interpreted as quasiparticle 
occupation numbers. With the help of the eq. (2.16) one can invert the relation (2.17) as 



(3.21) 



On the other hand, one can now express the field operator in terms of the natural orbitals 

as 

(/)(t, x) = {(f){t, x)) + J2 Kit, x) da{t) + <(t, x) dl{t)] (3.22) 

a 

with 

Vait.x) = Y.[fk{x)Uua{t) - .fl{x)VUt)] . (3.23) 

k 

Equations have been used to obtain these results. Therefore, our treatment for the general 
non-uniform field configuration involves expanding the field operator in the general natural 
orbitals, which are in turn given in terms of plane wave (or some other convenient) expansion. 

Ill-b. Formal Equations of Motion for the Gaussian Observables 

The next step is to derive the equation of motion for these simple variables {A^ , Pa , 
Uka and Vka)- For Ak{t) one finds immediately from the Heisenberg equation of motion 

iAk{t) = 1i[ak{t),H]F , (3.24) 

where H is the field Hamiltonian. The equation of motion for the remaining quantities 
can be obtained by taking the time- derivatives of the eigenvalue equation (2.12), again in 
close analogy with (2.13)-(2.14). Using the definitions for the matrices iV, X and i?, 
this equation can be written explicitly as 

^ WAU + WtlV + V^t[*U + V^A*V t/UV'* + WtlU* + l/tfiy* + V^A*U* ^ 



\ WA*V + U^tlU + V^WV + V^AU U^AU* + U^WV* + V^YiU* + V^AV* ) 

' P + [f/tf/ - VW , P\ ytf/* _ u\y* + {ytf/* _ f/ty* ^ p}^ ^ 

1 

^V^u -vW ^{v^ij -u^v ,p}^ p^\[j'^ij* -vW\p\ J 

(3.25) 
12 



where { }+ indicates an anticommutator. Note that there are two independent matrix 
equations, the remaining two being their complex conjugates. The block matrix elements 
in the left hand side of eq. (2.27) can be rewritten using the Heisenberg equation of motion 

as 



u^ku + u^tiv + v^ii*u + v^tv 



al3 



Tr 



(X/3 (X(y 5 ±1 



(3.26) 



U^kV* + U^tlU* + V^tlV* + V^k*U*] = Tr [dp d^ , H] F 



a/3 



(3.27) 



Equating corresponding block matrix elements on the two sides of equation (2.27) yields 



[P+[U^U -V^V ,P\] = Ti [d^d^^H] F , 



(3.28) 



v^u* - u^v* + {v^if* - u^v* , p} 



a/3 



Tr [df,d^,H] F . (3.29) 



Since P is a diagonal matrix, equation (3.28) splitts into two equations, one for a = f3 
and other for a j^ P , 



iPo 



Ti[did^,H]F , 



(3.30) 



t{P/3-Pa){U^U-V^V) 



afi 



Tr 



R ^ct 1 -^ 



aT^(3 



Moreover, equation (2.31) can be rewritten as 

I iy^W + U^V*) (l+Pa+P/3) = Tr [dpd^^H] F . 



(3.31) 



(3.32) 



For the particular case of a spatially uniform system, eq. (2.33) is trivial since the plane 
waves are the natural orbitals. 

In order to treat the dynamics of natural orbitals, which is given in terms of the awkward 
looking combinations UW—VW and VW* — UW* in eqs. (2.33) and (2.34), it is convenient 
to first define matrices h and g as 



h 



U^U - V^V 



(3.33) 



g = i (v^il* - U^V*) 



(3.34) 
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The next step is to find simple relations between U , V and h, g . The crucial point is to 
observe tat the dynamics of the eigenvector X can be described as being generated by a 
dynamical matrix Q such that 

iX{t) = X{t)n{t) G . (3.35) 

Using (2.16) and (2.17) one can solve eq. (2.37) for Q(t) as 

n{t) = iGX\t)GX{t)G = \ ^ • (3-36) 

\g* h* J 

Equating the matrix blocks of the equation (2.37) one finds 

iU = Uh + V*g* (3.37) 

iV = Vh + U*g* . (3.38) 

Therefore, the final equations of motion for the simple variables of the field are eqs. (3.24), 
(3.30), (3.37) and (3.38). The ingredients of the matrices h and g of eqs. (3.37) and (3.38) 
are found from (3.33) and (3.34). These equations are of course not closed equations since 
they still involve the full time evolution of the field operator as we have seen in previous 
section. 

The basic idea of our approximation scheme is described in section III. The additional 
difficult because of the spatial dependence of field can be handled when orbital representation 
is used. The crucial point here is to notice that the reduced density F(0)(t) in mean-field 
approximation can be conveniently written as 



Fo{t) = n 



l+Pa{t) 



Pa{t) 
l+Pa{t) 



dUt)da{t) 



(3.39) 



From these ingredients one can construct the correlated density F'{t) and finally the equa- 
tions of motion for a specific model. 

IV. Mean-Field and Collisional Dynamics in 0"^ Field Theory 

Thus far we have presented a general procedure to investigate kinetics of one-body observ- 
able in the context of a scalar field, without mentioning, however, any specific field model. 
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In this sections we will illustrate a specific example as an application of former more formal 
treatment. We consider the (h^ Hamiltonian 



H 



dx Ti 



H 



TT" 1 

y + 2 



^■^+i.(f%'l>?+'-^<l''+T.'l>* + 



m 

y 



9_ 

4! 



6m'^ 



(4.1) 
(4.2) 



The renormalization of (p^ theory in 1 + 1 dimension is well known ||T3| and can be achieved 
by introducing a mass counterterm drn^ given as 



Sm 



_9_ 
AL 



E 



1 



yW^+rn? 



(4.3) 



The expansion of the hamiltonian on the basis of natural orbitals follows directly from the 
discussion of section II. Therefore we have now all the necessary ingredients to implement the 
proposed approximation to the collisional dynamics. This is a lengthly by straightforward 
algebric exercise. The resulting equations of motion are 






|(^^fci+A:2+fc3+A:,0 ^l^ ^fcj ^L + ^fci+fea+fcs-fe , ^fci ^fca ^fcg^ 
3 (^4i+fc2-fe3+fc,0 ^fci ^li ^H + <^fci+fc2-fc3-A:,0 Ai M2 ^fcajl 



8L 



k\k2kz 



V ^01 ^02 ^03 ^0 



(•^fci+fe-fcs.O ^fei + 4i-A:-fc2+fc3,0 ^fcij zZ '^k2aTl^^{l + 2pa 



8Lko 



{A, + A*_, J2 



Pa = ^p{t) , 



+ ^^A(t) 



(4.4) 



[l+Pa+Pp)g, 



+ 



E 

k L 

_9_ 
8L 



aP 



k"^ + m? 
2ko 



V^ka ~ ^-ka)W-kl3 " ^kg) 



E 



kik2k3k4, 



l+Pa+P(3 
V ^01 ko2 ^03 ko4 



|(4i+fe2+fc3+fe4,0 ^fcl ^fc2 + ^ki+k2-k3-k^,0 Ak^^ki 
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9 



+ 



8L 



2^ 4i+A:2+fc3-fc4,0 / ,^ ,^ ,^ ,^ Tl^aTk2l3 2^ ^4 «' -^^4 a' ( ^ + 2Po 



kik2k-j.k4 



9 V- 1 + Pa + P/3 ^* ^* V- 



V ^01 ^02 ^03 ^04 
1 



8L 



h 



01 



fc2 



/co + m'^ 



zTg{t) 



(4.5) 



(P/3 - Pa) Kp 



+ 



E 

9 



16L 



■ E ^ 

feifc2fe3fe4 V^Ol ^02 /i^OS ^04 



{PP - Pa 



+ 4i+fc2-fc3+fc4 ,0 (^fci ^fca ^Xga Tk^fS + Al_^ Al^ TksfB T^^a) 

+ 26k^-k2+k3-k4,o [A*kiAk2Tk^aTk4/3 + AkiAl^Tk3f3T^^a)\ {Pp-Po 

9 V^ '^fei+A:2-A:3-A:4,0 



8-^ fcifc2fc3fc4 V^Ol ^02 hz ^04 



^ 



1 



8^ fci ^01 '" ' fc2 \/A;i + m2 



^fcl Q 2"fc3 /3 ^ Tfc*2 „/ Tfc^ q/ (P/3 - Pa 
a' 

1 



(P/J-Pq) +^rfe(t) , 



(4.6) 



where 



T, 



ka 



Ukr, — Vi 



ka 



and the colhsion integrals r(t) are 



^A{t) 



9^ 



96L2 



E 

k\k2k-i 



ykoi ^02 A;o3 A;o q^ 



E 



Q2Q!3 



'-'fci+A;2+fc3+fc,0 -'fciai -'fc2a2 fc3Q3 -'010203 

X 1^ T^ 1^ ji'^)* 

^ki+k2+k3-k,0 -ffcioi -^^202 -'-kia-i -'010203 
+ <JOfc^-|-fe2+A:-A:3,0 -'feioi -'fc202 -'^303 -'010203 
~ <jOfci+A;2-A;-fc3,0 -t/cioi -'fc202 -'^303 -'01O203 j 



(4.7) 



rp(t) 



96L2 



E 

k-ik2k-j.ki 



1 



V ^01 ^02 ^03 ^04 



E 

010203 



r ry-i* rri* rri* rp^ T\^) 

0A;i+fc2+fc3+fe4 ,0 -'feioi -'A;202 -'^303 -' fe40 -'0102030 
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I or rj-i^ ry-i* rji^ rp /'(2) 

"T O0ki-^k2+kz-k4^ , J- kiOLi -^ k2CX2 ksa -^ k^a^ ^aic 



a2CtOi3 



r ATI* ry-i* ry-1* /tt r(2) 

~ ^/i;i+fc2+/i;3-^4 , ^ kiai ^ k20L2 k^az ^4^ ^cyiOi20izo^ 
+ ^f^/i;i+/^'2-/c3-/c4 , -^ fcia -^ ^2^1 ^3^2 -^ /i^4a3 -'aaia203 J 

32-f'^ fcife2fe3fc4 V^Ol ^02 A;o3 A;o4 aiaj 

1 (^'^fei+fc2+A:3+fc4 ,0 ^fci + '^-A:i+fe2+A:3+fc4 ,0 ^fci j ^fcaoi Tk.>.a2 Tk^a ^aia2a 

2 (^'5fci-fc2-A:3+A:4 ,0 ^fei + ^ki+k2+k3-k4 ,0 ^li) ^feai -^fcsa ^^4^2 ^aiaa2 
— (4i-fc2-fc3+fc4 ,0 ^fei + ^ki+k2+ka-k4,oA*ki) '^k2ai ^4^2 -^^40 ^aia2a/ 

+ c.c. , (4i 

^2 1 

^^^^ kik2kik4 V^Ol ^02 ^03 ko4 aia2a;i 



0fci+A;2+fc3+fc4 , -f fciai -^ fc202 fc3a3 -^ fc4Q -"q 



r{l) 
i+A;2+fc3+fc4 , -L fciai ^ fc202 "^ ^303 "^ k^a ^aia2a;if3 

+ <J"fci+A;2+A;3-A;4 , -^ fc^Qi -^ ^202 ^3^ ^=403 -'aia2/3Q3 

r 7^ 7"" 7^ 7""* r^^)* 

t'fci+fc2+fc3-fc4 , -t feiQi J- k2CX2 -L k^a^ J- k^a -'010203/3 

I Q 'T ry-T* rp^ rp rp T\^) 1 

+ ^Oki+k2-k3-k4,,0 -'fciQ-'fc2ai -'fc3"2 -'- k^cii J- paia2a3 J 

32-f'^ fcife2fe3fc4 V^Ol ^02 ^03 ^04 aia2 

|('^-fci+fe2+fe3+fe4,0^fcl + '^fci+fc2+A:3+fc4,0^fcij ^feoi ^fc3a2 -^^4" -^aia2/3 

2 ( (5^1-^2-^3+^4 , ^fei + '5fci+fc2+fc3-A:4 , ^fci ) ^fc2ai "^ksa Tk4a2 ^ail3a2 



— [^5ki-k2-k3+k4,oAki+Ski+k2+k3-k4,oA*ki) ^feai ^fc3a2 ^fctia -^aia2/3/ 

+ C.C. (a^/3) , (4.9) 

Yg(t) = ^ J2 ^ E 

^^^^ kik2k3k4 V^Ol ko2 A;o3 A;o4 aia2a3 

/a 1^ T T 7~'* r^^)* 

|"fci+fc2+fc3-fe4,0 -tfciai -tA;202 -'A;3a3 -^ fc4Q -'aia2a3/3 
+ <J"fci+fc2-A:3-fc4 , -f fcioi -^ ^202 ^ k-iOt ^ k4a:i -'aia2/3Q3 

r rj-i* ry-i* ryn* ^* T\^) 

~ 0ki+k2+k3+k4 , -L ^jQ,^ J ^2a2 ^303 ^40 -'aiQ2Q3/3 
"T dyfe-^+fc2+A;3-A:4 , -^ fciQi -^ ^2^2 ^30 ^4^3 -'aia2/3Q3 J 

+ ^ E ^ E 

22-^^ fcife2fe3fc4 V^Ol ^02 ^03 ko4 aiaa 
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1 [0ki+k2+k-i-k4,0 ^ki + 0ki-k2~k3+k4 ,0 ^ki) ^feai ^fc3a2 ^A:4a ^aia2P 

2 {6k-^-k2-ki+k4 ,0 ^ki + '^fci+fc2+fc3-A:4 ,0 ^fci ) ^fcaai ^fcaa ^fc4a2 -^ai/3a2 



5-fci+fc2+A:3+A:4 ,0 ^fci + <^fei+fc2+fc3+fc4 ,0 ^fej ^fc2ai ^4^2 -^^40 -^aia2/3| 
a ^ P . (4.10) 



In these equations we used the abbreviations 

{-i\ If X — ^ ^k^ 4-k' A-k! 4-k' x — ^ / 

•"" fc^fc^fe^fc^ \l Hi 1^02 Hz l^m 71727374 

4 4 4 

1 + II P7, + H P7. P7, + H Plr Pi J Pit 

i i<j i<j<(. 

(m:^^^ (t, t') m:^^^ {t, t') m:3^3 (t, t') m:^^^ (t, to) , (4.1 1) 

^Sa2a3"4W = /, ^^' H /// ,'/ w "r r ^ (^^•'i7i^fc^72^fc^73 ^474)4, 

'^ k[k'^k'^k'^ \/HlH2H3H4 71727374 



P74 P71 P72 1^73 + P74 2^ Pli + P74 Z^ Pli Pli 

(M:^,,{t,t') M:^^^{t,t') Ml^^.^{t,t') M^,,,it,t') ) , (4.12) 



I^aLa^aAt) - Jdt'j: llTlX^^U, „, ^ „ (^^•'171^^^72^473^474),, 



^k[+k'^-k[^-k'^,0 
'° k[k'^k'^k'^ yKl ko2 ko3 Ka 71727374 

[1)73^74(1 +P7i)(l+P72) - ^71^72(1+^73) (1+^74)]*' 

{m:^^^ (t, t') M:^^^ {t, t') M«3,3 (t, t') M,,,, (t, t') ) , (4.13) 

^Sl2a3W = dt' Y. 1., ., ., = E (^fc>^fc^72^fel73),, 

•^° k[k'^k'^k'^ y koi ko2 kQ3 kQ4 717273 



1 + II P7. + H P7. P7, 

(M*^,^(t,t')M:^^^(t,t')M:3^3(t,t')) , (4.14) 



•'^ k[k'^k'^k'^ \jKlK2KzKA 717273 

\^k[+k'^+k'^-k'^,Q M[ + h[-k'^-k'^+k'^,0 ^k'J ^, 
[1^73(1+^71+^72) - ^71^72]*' 

(M:^^^(t,t')M:^^^(t,t')M„3,3(t,t')) . (4.15) 

In summary, we have presented in this paper a framework to study real-time evolution 
of scalar field theory. The technique has been applied to nuclear many-body theory and 
extended recently in the context of homogeneous field configurations. Here we discuss this 
general problem when the spatial dependence is important. Here we show that the spatial 
dependence can be treated in the general orbital representation and the dynamics of these 
orbitals are expressed in a closed form by a set of selfconsistent equations. In this way, 
the time-dependent technique can be used to improve the usual gaussian-like mean field 
approximation, where the coUisional dynamics are given by approprite memorial integrals. 
We have illustrated these procedures within the simplest context of self-interacting 0^ theory. 

Appendix A: Projection Technique in Field Theory 

In order to calculate equations of motion (2.24), (2.32), (2.33) and (2.34) we first decom- 
pose F in two parts 

F = Foit) + F'it) (A.l) 

where i^o(^) is the exponential of a one-boson density given in (3.39). A crucial point is to 
observe that -fo(t) can be seen as a time-dependent projection of F , i.e.. For the explicit 
construction of P(t) we require, in addition to eqs. (3.3), the condition 

iFoit) = [Fo{t),H] +]P{t)[H,F] (A.2) 

which is the Heisenberg picture counterpart of the Schrodinger picture condition used in 
ref[pn| to determine P uniquely. The resulting form for P(t) is (see ref. for details of 
the derivation) 



P- 



did. 



^ai ^012 Pa2 ^aic 



y "«"» P» T^(.) + y ""i""^ P^2^a^a2 ^ /^t ^ 



Pa2(l+Pai) 
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E 



'"T.(4-) + ^T.w 



p. 



+ 






2(l+p«J(l+p,J 



Pa 

ir \U1a2 Q!i 



<ia-) 


+ >; 

«1«2 


) 


[ Fo . 



^^^^^ Tr (4^ 4^ 



(A.3) 



The next step is to obtain a differential equation of F'{t). This follows immediately 
from eqs. (3.1) and (3.4), 



\j^ -P(t)lL'j F\t) = Q(t)lLFo(f) 



(A.4) 



where we introduced the operators 



Q(t) = IL-P(t) , I.- = [H, 



This equation has the formal solution 



where G{t, t') is the time-ordered Green's Function 



(A.5) 



F'it) = G{t,{])F'{{]) - i t dt'm,t')<i^{t')l.Fo{t') , (A.6) 

Jo 



G(t,t') = T expi [ rfrP(T)L 
Jt' 



(A.7) 



A systematic expansion to treat the memory integral of the equation (3.8) has been 



discussed in ref.|Tl[ in the Schrodinger picture. The implementation of the corresponding 
expansion in the Heisenberg picture consists in approximating the time evolution of the field 
operators by the simpler mean-field Hamiltonian 



Ho = P^H + J2diTi[da,H]F'-Y.daTT[di,H]F' 

a a 



o-ia.2 "V-^ "I" Poii ~r 'P0L2) 

~ hi. 2(1 + p. 



A 0-2 



a\a2 



+ Pai + P, 



Tr 



Ot2j 



d^ d^ H 

"ai ""02 ' 



F'. 



(A. 
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Using this approximation, one can solve explicitly the Heisenberg field operator equation as 

d^{t) ~ Y. M^^{t,t') d^{t') , (A.9) 

7 

where the matrix Ma^ is the solution of the matrix equation 

iM{t,t') = F{t)M{t,t') , (A.IO) 

and the matrix F involves matrix h and mean-field energy (see appendix A for details of 
the derivation). 



In this way, the authors of the ref.|ll| devised a systematic expansion for the correlation 



density F'(t). The corresponding expansion in the Heisenberg picture is 

F'{t) = Git,0)F'iO) - i f'dt[Q{ti)LFo{ti) 

Jo 



dti 





dt2Q{h){L-LQ{t2) 



Q(ti)LFo(ti) + --- , (A.li: 



where Lq ■ = [Hq , •]. In what follows we restrict ourselves to initial conditions such that 
F'(0) = and to the lowest approximation for F'{t). Therefore, the evaluation of the 
equations of motion involves traces of the type 

TT[d{t),H]Foit) - 2Tr[6(t),i/] f dt' qit') [H, Fo{t')] , (A.12) 

Jo 

where 0{t) can be da{t) , d'^^it) dp{t) , daif) dj^it) , and operators at different times are 
related by equation (3.11). 

Appendix B: Approximation for the Time Evolution of 
the Heisenberg Field Operator 

In section II, we have discussed that our approximation consists in replacing the time 
evolution of the field operator by a simpler mean-field Hamiltonian given by equation (3.10). 
We show now that this allows one to solve the Heisenberg operator equation 

id^ = [da,Ho]-J2 f^»ida-J2 9-/ad\-Tr[da,H] . (B.l) 



7 
7 
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The last three terms are the exphcit time dependence of da{t) related to shifted amplitudes 
Ak{t) and to effects of the general Bogoliubov transformation for operators (2.19). 

First, we write P^H , using ciclyc properties of the traces, as (see ref. [ ] for the 

construction of P^) 

P'H = fl-E ^^^:^^] Tr (H Fo) + ^ < ^^-^ " ^"^ ^^ Tr (4^ d., // Fo) 

- Ed„Tr[4,if]Fo + 5]4TrK,iJ]Fo 

a a 



ir'^ 2(1 + p. 



~r Pai + Pa2j 



dl dl 



i:t„ 2(1 + p. 



02 



+ Pai + Pa2j 



Tr [4^ 4^ , H] Fo 



Tr [d«2 rfc^i , H] Fq 



(B.2) 



Hence, equation (3.10) becomes 



Hn 



dl. dn 



p 



1-E -;;/" Tr(//Fo)+ E 



ai«2 



Pa2(l+Pai) 



Tr (rfL da. H F 



"0.2 "^«1 



■ai n ro) 



- Y.daTT[dlH]F + J2diTT[da,H]F 

a a 

_ ^ da,da2 ^ Tr[4^4^,i7]F 

aia2 ■^{^ ' Pai ' Pa2) 

d^ d'l" 

_|_ y^ "ai "^2 

«i«2 ^ I -'- "T Pai "T Pa2 y 



Tr [da2 (io,^ , H] F 



(B.3) 



Using (A. 3) in (A.l) one obtains immediately 



idn 



dn' 






(B.4) 



The calculation for the first two terms is straightforward. It yields 



a' k 



ko {Uka + V*ka){Uka' + V^ka') + 7 ^fca 7),^.' 



(in 
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+ JL- y: e 

^6^ kik2k3k4 V^Ol ^02 ^03 ^04 (,/ 

X |4i+fc2-fc3-fe4,0 (A'l ^fc2 ^4" ^^40' + ^fcl ^L ^fcsa' ^fc4Qy 
+ 5fci+fe2+A:3-A:4 ,0 (^fci ^fe ^fcsa' Tj^^a + ^fci ^fca -^^aa ^^40'^ 

+ 24i-fc2+fc3-fe4,0 (^fci ^fc2 ^4" ^fc4a' + ^fci ^fc2 ^fcsa' ^fc4a j I ^"' 

_|_ ^ V^ ^kl+k2-k3-k4 ,0 ■^^ ry-,* ry-,* 1^ T rl 

^ AT 2^ ft. T, T, T— 2^ -^ fcl a -^ fc2 02 -^ ^3 ai -^ fc4 "2 Pq2 "qi 

^^ kikik-iki V«^0lfc02fc03fc04 aiQ2 

_|_ 9 Sr^ Oki+k2-ka-k4, ,0 ■^ rp* rp* rp rp rl rl ^ 

+ or Z^ /T r r TT- Z^ -^ fci ai -^ A:2 02 -^ '=3 "2 -^ A:4 ai "ai "qi . 

^-^ fcifc2fc3fc4 V«^0lfc02fc03/C04 Q1Q2 

^^ kik2 Jkl + m? 1^02 a. 



= E 5«7^7 • (B-5) 

7 

One finds finally the Heisenberg operator equation as being 



I 



(Xfj — / -LJ (y^ry LLfy 7 ' '-'(yy ^' 



7 



7 7 



E ^"7^7 (B-6) 



7 



or in the matrix form 

id{t) = F{t)d{t) . (B.7) 

It is easy to see that the solution for (A. 7) is 

d{t) = M{t,t')d{t') , (B.8) 

where M{t, t') is the solution of the matrix equation 

tM{t,t') = F{t)M{t,t') (B.9) 
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with initial condition 

M{t,t') = 1 . (B.IO) 
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